Color singlet suppression of quark-gluon plasma formation 



Jes Madsen, Dan M. Jensen and Michael B. Christiansen 

Institute of Physics and Astronomy, University of Aarhus, DK-8000 Arhus C, Denmark 

(June 7, 1995) 



G^ . Abstract 

' The rate of quark-gkion plasma droplet nucleation in superheated hadronic 

^ ' matter is calculated within the MIT bag model. The requirements of color sin- 
gletness and (to less extent) fixed momentum suppress the nucleation rate by 

' many orders of magnitude, making thermal nucleation of quark-gluon plasma 

o ■ 

\ droplets unlikely in ultrarelativistic heavy-ion collisions if the transition is 

^ \ first order and reasonably described by the bag model. 
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Ultrarelativistic heavy-ion collision experiments at CERN and Brookhaven aim at the 
formation of quark-gluon plasma. No unambiguous signals have been detected so far, but 
the prospects will improve tremendously with the next generation of colliders (LHC and 
RHIC). Normally it is assumed that formation of a quark-gluon plasma will take place if 
a region in the colliding nuclei reach a combination of temperature and chemical potential 
that brings it into that part of the phase diagram for hadronic matter where the quark-gluon 
plasma has a lower free energy than a hadron gas. 

Should the phase transition be of first order, the rate for quark-gluon droplet formation, 
TZ, can be estimated in the framework of homogeneous nucleation theory , where 

7^ ^ T^exp(-AF/r) ^ 0.2im-'\lO~^hy^T^oo(^M-^F'/T). (1) 

Here AF is the height of the free energy barrier which the thermal nucleation has to over- 
come, Tioo is the temperature in units of lOOMeV, and the rate has been expressed in units 
of typical heavy ion collision volumes and time scales. There is an extensive literature on 
the proper choice of prefactor, but since the exponential is by far the most important for 
our present discussion (the prefactor does not influence the relative rate suppression due to 
the effects discussed in our paper), we have chosen the dimensional estimate T^. We note 
that a more realistic choice of prefactor |0] could change the nucleation rate by a few orders 
of magnitude 0. A quantitative estimate of the effect in the present context is hindered 
by the fact that no derivation of the prefactor exists for a situation where curvature rather 
than surface tension is the most important contribution to the surface energy. 

It is normally assumed that a moderate superheating is sufficient to allow droplet nucle- 
ation. This is confirmed by calculations of the droplet free energy in simple phenomenological 
models like, e.g., the MIT bag model Such calculations have, however, neglected the 
requirements of color singletness and fixed momentum for the quark droplet. Both of these 
constraints significantly reduce the effective degrees of freedom in the quark-gluon plasma 

thereby increasing the free energy barrier, AF, and reducing the nucleation rate. As 
we show below, the nucleation rate is in fact reduced by many orders of magnitude, making 
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thermal nucleation virtually impossible regardless of the amount of superheating. 

We note that the ingredients of our rather simple calculation have been around in the 
literature for about a decade, but to our knowledge no-one has checked these dramatic 
consequences explicitly before. 

We want to calculate the nucleation rate self-consistently within the MIT bag model, 
including two flavors of massless quarks, gluons, a bag pressure, volume and curvature 
energies, and the constraint that quark droplets must have a fixed momentum and be color 
neutral. It will turn out that our conclusions are not very sensitive to chemical potential, 
bag constant, number of quark flavors or the inclusion of the hadron gas (in most of the 
following we therefore just ascribe a chemical potential but no pressure to the hadrons), but 
very sensitive to the color singlet constraint. 

Elze and Greiner have derived the color singlet fixed-momentum grand canonical 
partition function, Z(T, i?, /x; p), for a hot, spherical quark-gluon plasma droplet of temper- 
ature T, radius R, quark chemical potential fi (corresponding to a baryon chemical potential 
fiB = 3/i), and total momentum p. From the partition function one can find the grand po- 
tential, Q = —TlnZ, as 
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Nq is the number of massless quark flavors (we take Nq = 2 in most of the following). V is 
the volume of the droplet. The first term in the grand potential is the contribution from the 
color singlet projection. The next two terms are from fixing the momentum (notice that the 
first of these terms survive even for p = 0), and the last two terms are, respectively, the bag 
energy proportional to the bag constant, B, and the "normal" bag model grand potential 
of the quark-gluon gas. Notice that finite size effects, i.e. energy terms proportional to the 
extrinsic curvature (Svri?) of the droplet have been included, whereas massless quarks and 
gluons contribute no surface tension term proportional to the bag area. One could add to 
Eq. (1^) a surface energy from the hadron phase, as well as the hadron pressure times volume 
(i7 1-^ i7 + Pnadron^ + ^Surfaco.Hadron, whcrc Pnadron IS the prcssurc of the hadrou phase). 
These terms would both add positive contributions to Q, thereby further decreasing the 
quark-gluon droplet formation rate. For reasonable choices of parameters these terms are 
not decisive, and since our aim is to clearly demonstrate the enormous rate suppression due 
to the color singlet and momentum constraints, we have chosen not to incorporate them. 
Had we done that, they would only further suppress the nucleation rate and strengthen our 
conclusion (we return to this issue below). 

The partition function was derived in a saddle-point approximation, which is expected 
to break down for RT 0. When /x = it should be good to 30% for RT = 1 and a few 
percent for RT ^ 2, and the error decreases rapidly for /i > 

The grand potential as a function of droplet radius for fixed bag constant, fi = 0, and 
a temperature corresponding to nearly 30MeV superheating is shown in Figure 1. Curves 
are given with and without inclusion of the color singlet and p = constraints. One notes 
a significant increase in Q, in particular due to the color singlet requirement. 

In Figure 2 we show the corresponding nucleation rate of quark-gluon droplets estimated 
according to Eq. (|1|). The change in free energy, AF, necessary for formation of a critical 
bubble in chemical equilibrium with the surrounding hadrons at fixed temperature is just the 
height of the Q barrier, so AF = Q with Q given in Eq. (^. One recognizes the well-known 
reduction of the transition temperature for increased chemical potential, but what is more 
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important in the present context is the significant suppression (by 4-5 orders of magnitude 
at extreme superheating, and much more at moderate superheating) of the nucleation rate. 
Units in the plot are chosen such that unity corresponds to 1 nucleation per fm'^ per 10~^'^s 
(=3 fm/c). Typical volumes and time scales relevant for quark-gluon plasma formation in 
ultrarelativistic heavy-ion collisions could amount to maybe 100 fm^-fm/c, so in standard 
calculations one predicts a fair probability for bubble nucleation in the case of, say, a 
lOMeV superheating. However, the suppression resulting from requiring fixed momentum 
and color singletness gives a very low probability, even for (unrealistic?) super heatings of 
50-100 MeV. 

This surprising conclusion is not significantly influenced by changes in the assumed num- 
ber of massless quarks (here taken to be 2, corresponding to infinite s-quark mass), or to 
the external pressure contributed by the hadron gas (the latter corresponds effectively to 
a small increase in B 0). Other choices of B would correspond to a rescaling of T and n 
in Fig. 2 proportional to B^^^, and of TZ proportional to B (i.e., no qualitative change in 
the conclusion). A massive s-quark would add a surface tension term to the energy, and a 
surface tension contribution could come from the hadron phase as well. This would change 
the numbers (always in the direction of even lower nucleation probability and a slightly 
higher critical radius), but not the effect of color singlet suppression. No derivation exists 
of the color singlet partition function for a massive quark; neither does a calculation for 
non-zero strong coupling constant, which has therefore been assumed equal to zero in the 
present investigation. 

The results do of course depend on the choice of model for the quark-gluon plasma. 
The MIT bag is certainly relevant only if a first order phase transition is involved at all. 
Some lattice calculations indicate, that the transition could be second order for yU = 0, but 
the issue is far from settled, and calculations for /i 7^ do not exist. Other models would 
give other numerical suppression factors, but the color singlet constraint will in any case 
reduce the effective number of degrees of freedom thereby increasing AF in any model, 
so qualitatively a rate suppression should be expected. 
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Equation (0) is based on a saddle-point approximation, which is again based on an 
expansion of the density of states in terms of volume and curvature terms. The expansion 
is known to reproduce direct summation very well for yU = [|^ and for T = 0], so there is 
no reason to disbelieve this assumption. The uncertainty in the saddle-point approximation 
was discussed above. It is not negligible (but not devastating either) for critical bubbles with 
RT ^ 0.9 as typically found for superheating beyond 20MeV when /i = 0, but negligible for 
the much larger bubbles involved closer to the bulk phase transition temperature, and for 
fi ^ where the saddle-point approximation is much better p. Figure 3 shows the radius 
of critical bubbles as a function of temperature. 

Thus, in spite of all the reservations, we conclude that a very significant suppression of the 
nucleation rate for quark-gluon plasma droplets in a superheated hadron gas is an inevitable 
consequence of the fixed momentum and color singlet constraints if the quark-hadron phase 
transition is first order. Other mechanisms than thermal nucleation (e.g., nucleation due 
to impurities) are needed to form quark-gluon plasma droplets in ultrarelativistic heavy-ion 
collisions if the phase transition is first order. 
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1 + 0.33(T/S^/4)4 
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FIGURES 

FIG. 1. The grand potential, J7, as a function of droplet radius for S^/^ = 200MeV, /x = 0, 
and T = 170MeV. The lower curve is without momentum and color singlet constraints, middle 
curve with the color singlet constraint, and upper curve with both constraints included (at zero 
momentum). 

FIG. 2. The nucleation rate of quark-gluon droplets as a function of temperature for 
^1/4 = 200MeV. Sohd curves for quark chemical potential zero, dotted curves for n = lOOMeV, 
dashed curves for = SOOMeV, dot-dash curves for /x = 400MeV, and dash-triple dot curves for 
H = 500MeV. The lower set of curves includes the color and momentum constraints; the upper set 
does not. 

FIG. 3. Radius of critical bubbles as a function of temperature for calculations including color 
and momentum constraints, for chemical potentials as in Fig. 2. 
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